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We evaluate the levels of the rnr system in the p channel in finite volume. We investigate the 
dependence on the ir mass and compare with other approaches which use QCD lattice calculations 
and effective theories. We also illustrate the errors induced by using the conventional Liischer 
approach instead of a more accurate one recently developed that takes into account exactly the 
relativistic two meson propagators. Finally we make use of this latter approach to solve the inverse 
problem, getting nn phase shifts from "synthetic" lattice data and show which accuracy is needed 
£Nj , in these data to obtain the p properties with a desired accuracy. 
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I. INTRODUCTION 



The determination of hadron spectra is one of the challenging tasks of Lattice QCD and many efforts are being 
devoted to this problem oJ— FlQll . As one is aiming at higher accuracy and consistency with the decay channels of the 
hadrons, Liischer's approach [20l[2l|. reconstructing phase shifts of the decay products from the discrete energy levels 
of the box, is playing gradually a more important role. Other strategies formerly used, like making use of the "avoided 
O j level crossing", are leaving room to other methods once this criteria has been shown insufficient for resonances with 
<D . a large width [^Ull. 

One of the hadrons that has received more attention is the p meson. There have been a few attempts to describe 
the p using Liischer's approach [25U30l ]. On the other hand the first lattice estimate of the p — > tttt amplitude [3l[ did 
not apply Liischer's method. The most recent work on the p properties from the lattice point of view is the one of 
k*" ! and we shall refer to it for comparison along the paper. In between in [24| Liischer's approach has been recently 
simplified and improved by keeping the full relativistic two body propagator (Liischer's approach keeps the imaginary 
part of this propagator exactly but makes approximations on the real part) and extending the method to two or more 
coupled channels. The method has also been applied in to obtain finite volume results from the Jiilich model 
for meson baryon interaction and in [33$ . to study the interaction of the DK and r]D s system where the D s *o(2317) 
resonance is dynamically generated from the interaction of these particles [3~4T - l37j . The case of the k resonance in the 
Kir channel is also addressed along the lines of [23 in (38| and an extension of the approach of [24[ to the case of 
interaction of unstable particles has been made in [39 L 



In the present work we shall use the approach of 2J| to the case of the interaction of two pions in p-waves and 
isospin 7=1, the p channel. On the first part we shall use the chiral unitary approach to 7T7t scattering in p-waves 
from [4(| 5l| and apply it to get energy levels in a finite box. We shall then see how they change with the pion mass 
to establish connection with lattice results which run for heavier masses than the physical pion mass. In a second 
step we will face the inverse problem of getting tttt phase shifts from the energy spectra in the box, determining the 
precision needed in the lattice spectra to get the phase shifts with a demanded accuracy. 



II. FORMALISM OF THE CHIRAL UNITARY APPROACH FOR THE p IN INFINITE AND FINITE 

VOLUME. 

We shall follow here the approach of [4l| to get the 7T7T scattering amplitude in p-waves. In this work the pion 
and kaon form factors were studied using the chiral Lagrangians of [42| and those for the coupling of resonances to 
pairs of pseudoscalar mesons of [43} . The work was done obtaining a kernel, or 7T7t potential, from those Lagrangians, 
which contains a contact term from the lowest order Lagrangian of [42| and another term, the explicit p pole, with the 
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coupling of the p to the pions obtained from [43| . One includes also the KK system and unitarizes the two channel 
problem using the coupled channels Bethe Salpeter equations, or the equivalent N/D method. Although one can 
approach the problem from different perspectives, like the inverse amplitude method (IAM) [HI, S^|, or the Bethe 
Salpeter equations with an energy dependent potential provided by the first and second order terms of the chiral 
Lagrangians of [42|, as done in [4(31 ] . the approach of [4l| by separating the contact term and a bare p pole, already 
accepts from the beginning the proven fact that the p is not a tttt resonance, but rather a genuine state which has 
tttt as decay channel. This has been checked thoroughly by looking at the large N c behaviour of the pole position in 
[13,13, where the mass of the p stays constant as N c grows and the width decreases. This is opposite to the composite 
states, dynamically generated by the interaction, like the light scalar mesons, where the width grows considerably 
as N c increases. A different perspective is provided in [4^], where it is shown that all chiral logarithms cancel out 
in the rho-channel, while they do not cancel for the sigma case, and they strongly influence the properties of this 
latter resonance. In other words, the loops from intermediate tttt states are very relevant for the sigma and highly 
unimportant for the p case. In even simpler words we could say that the sigma is a tttt resonance, while the p is a non 
tttt state which decays in this channel. One way to take this into account is to allow for a CDD (Castillejo, Dalitz, 
Dyson) pole in the potential [50], and this is what naturally appears by combining the Lagrangians of 1421 and W3J 
and is implemented in the approach of [4lj . 

In the chiral unitary approach one uses the coupled channel Bethe-Salpeter equations in their on shell factorized 
form (Hi [Hi 



T=(1-VG)" 1 V. (1) 
The relevant ^-matrices for 7T7t and KK scattering have been studied in Ref. pi] ]: 



,/=i ^ 2p»/ , 2G 2 V 



V I — 1 rp £_7T / i 



V I=1 \P)T V2p*PK f 1G\ S \ 



7T7r,7r7r -^71^7? ,7T^7T 3 f 2 V J 2 A'l 2 

T/ 7=l _ rp rp _ PK fl I 2( ^V S 

KK.KK ~ 1 K+K~,K+K- - J-K+K-,K«K« ~ ~Jp\ l + ~JT M 2 _ ~ 

,1=0 rp . rp _ PK fo j 2( ^V s , 4Gy 



V KK,KK ~ T K+K-,K+K- + T K+K~,K°K ~ ( 3 + ~JT~ M 2 _ ~ + ~~pT 

We note that these definitions have a minus sign difference with respect to [4l| to follow the more standard notation 
for V and T. We work in the isospin base, where the states are defined as: 

|vT7r) /=1 = i|/T+/T~ - 7T~/T + ) , (3) 

\KK) t=1 = ^\K + K- - K°K°) , (4) 

\KK) I=0 = -=\K + K~ + K°K°) . (5) 

Note the extra l/\/2 in the normalization of the pions (unitary normalization) which is adopted to account for the 
identity of the particles in the intermediate states. Since in this paper we study the p meson which has isospin 1=1, 
there are only two channels to be considered, \tttt)i = i and \KK)i—\. For simplicity, we rewrite the ^-matrix as 
V^^ = V I=1 Voo = V I=1 - and Vv> = V I=1 - 

"11 7T7T , 7T7T ' "22 KK,KK * 12 TTTT^KK 

The parameters in Eq. ^ are also taken from Ref. [4l| , in particular: 

G v = 53MeV , 

/ = 87.4MeV . (6) 

But for the bare p mass M p of Eqs. © we choose the value 837.3 MeV, which is slightly larger than the value 829.8 
MeV used in Ref. 41]. Doing this, the phase shifts obtained are slightly better compared with the experimental 
data [53] , as shown in Fig. Q] in the next section. 

The G-function for the two-meson propagator having masses m and M is defined as 

r( 2, ■ f d 4 g 1 1 m 

{P ' J (2tt) 4 q 2 - m 2 + ie (p - q) 2 - M 2 + ie ' { > 
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In Ref. |4l| it has been calculated in the continuum space (infinite volume) using dimensional regularization, and the 
result is: 



Gi(s) 



6tt 2 V 



2 + di + ctj(s) log 



<Ji{s) + 1 



(8) 

1) or \KK) I=1 (i = 2), and 

tTi(s) = y/T— ^mf/s. The parameters di are also taken from Ref. [4l| which are fitted with the experimental data of 
7T7T P-wave scattering: 



16tt 2 V " ' ' "•v-z-o _ Xj 
where the subindex i refers to the corresponding two-meson state, either |7T7r)/=i (i 
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(9) 



with p = M p . Eqs. ^ are taken such as to guarantee a perfect matching with the results of chiral perturbation 
theory for the form factor of p54j . 

All the above formulae are for the study of the p meson in the infinite space. To study the p meson in the finite 
volume, we just need to change the G-function of dimensional regularization (Eq. ([5])) by the one which is defined in 
the finite box of side L. This was deduced in [33j and it is given by G(s), defined through: 



G(s) 



9<9n 



d 3 q 



•Ai) 



(10) 



2 Wl (q> 2 (g) P 2 - + u 2 {q)f 



(11) 



where u>i^{q) — y rn\ 2 + and the discrete momenta in the sum given by q = -^-n [n G Z ). In this way we have 

just changed the integration over momenta by a sum over the discrete values of the momenta allowed by the periodic 
conditions in the box. 

The real part of G(s) — G(s) is shown in Fig. Q] as a function of (? max - We note that we have fixed L to be 2.5 
m~ 1 and E to be 770 MeV. We find that it has a good convergence when one goes to large values of q m ax- However, 
as already done in [U [H| it is more practical and equally accurate to make an average of this quantity for smaller 
values of q m ax with the consequent economy in computation. 
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FIG. 1: The difference between G(s) and Re[GDit(s)] calculated using Eq. (|10[1 as a function of q n 
and E = 770 MeV. 



Here we choose L = 2.5 m„ 



III. THE ENERGY LEVELS IN THE CHIRAL UNITARY APPROACH 

Using the G-function in Eq. (fit))) , we calculate the energy levels in the finite volume for different values of L, by 
looking at the poles of the T matrix, which appear when the determinant of 1 — VG is zero: 

Det(l - VG) = 1 - VuGxis) - V 22 G 2 {s) + (V n V 22 - V 1 \)G 1 {s)G 2 {s) = 0. (12) 
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As we have said before, our procedure follows closely the method used in Refs. [HI, H3, 38, 55], so here we go directly 
to the discussion of the numerical results. The energy levels for 7T7t P-wave scattering are functions of the cubic box 
size L, as well as the pion mass m n . We shall study the dependence on these variables. The volume dependence is 
shown here and the pion mass dependence in the next section. 

In the left panel of Fig. [2] we show the energy levels as function of L for the coupled channels which are obtained 
after performing an average for different g m ax values between 1200 MeV and 2000 MeV. On the right panel we show 
the same ones separately for values of <7 max = 1300, 1500, 1700 and 1900 MeV. We find that the energy levels for 
different g max are almost the same. At least we can not differentiate them in Fig. [21 This surprising result is simply a 
consequence of the fact that the p meson is very weakly tied to the nn loops, accounted for by the G function, where 
the L dependence appears. We note that the higher curves seem to cross each other. This does not happen because 
of the ordering chosen for the levels and then produces the "avoided level crossing", but they are a consequence of 
the couple channel dynamics. To see it clearly, we show the results for single irir and KK channels in Fig. [3J 




L [mi] L [mi] 

FIG. 2: Energy levels as functions of the cubic box size L, derived from the coupled channels unitary approach of Ref. [4ll ] and 
using G from Eq. (imp . We perform an average for different g max values between 1200 MeV and 2000 MeV in the left figure, 
while show them separately in the right figure for different g max = 1300, 1500, 1700 and 1900 MeV, which are almost the same. 



In Fig. [3] we show the energy levels for the single channel, either the 7T7t channel or the KK one. We find that the 
first and second energy levels in the left figure, which are are obtained using only the 7T7r channel, are quite similar to 
the results for the coupled channels. However, the higher levels are a mixture of both the irir and the KK channels. 
This is quite reasonable since the KK channel does not contribute much in the low energy region. Since the first 
(lowest) energy level should be well calculated using the chiral unitary approach which is very accurate in this energy 
region, in the following discussions we shall mainly concentrate on it, and consequently some of our calculations will 
be done only in the 7T7t channel. 




FIG. 3: Energy levels as functions of the cubic box size L for the nn channel and the KK channels, shown in the left and right 
hand sides, respectively. 



From these energy levels, we can obtain their corresponding phase shifts. To do this we use only the irir channel 
and follow the procedures used in Ref. [24[ ■ Since the energy levels are obtained by solving the poles of Eq. (fT2")l , they 
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satisfy the following relation: 



1 - VG ° nc channcl ) Vn^VniE)- 1 - G tl (E,L)) = 0, 



which can be used to calculate the scattering matrix: 



T n (E,L) = {VuiE)- 1 - GniE))- 1 = (G n (E,L) - G U {E)Y 



(13) 



(14) 



Since this relation holds for the energy levels shown in Fig. Eland Fig. [31 L and E are not independent. Moreover when 
we only consider the first (lowest) energy level, we have Tu(E, L) = Tn(E, L(E)) = Tn(E). From this calculated T 
matrix one can evaluate the tttt P-wave phase shifts which are related to our T by 



-8ttE 



pn cot 5(p) - ipu 



(15) 



In Fig. 01 we show the tttt P-wave phase shifts 8, calculated through three different approaches. The solid curve 
is taken from Ref. [4l| using the chiral unitary approach. This result has a good consistence with the experimental 
data [53|. The dashed curve is the phase shift extracted from the first energy level of the left figure of Fig. [3] obtained 
using only the tttt channel. The results are the same as the solid curve. The dotted curve is the phase shift extracted 
from the first energy level of the left figure of Fig. [21 obtained with the coupled channels but using only 7T7t in the 
analysis. We find that there are small differences in the high energy region, which are caused by the KK channel. 




400 500 600 700 800 900 1000 1100 



FIG. 4: The solid curve is the tttt scattering P-wave phase shift calculated using the chiral unitary approach [4l[. The dashed 
curve is the phase shift extracted from the first energy level of the left figure of Fig. [21 obtained using only the tttt channel. The 
dotted curve is the phase shift extracted from the first energy level of the left figure of Fig. [2] obtained using coupled channels. 



Using the phase shifts 6(E) we can fit the physical quantities for the p meson, such as m p , g p7r7r and T p . We note 
that m p is the p mass we obtained, i.e., one of our outputs; while M p is the bare p mass, i.e., one of our inputs. We 
use the following two equations to extract the p properties: 



cot 6(s) 



VsT p (s) 



and r p (s) 



P dflTT 



We note that the factor 8ir in the second equation is our normalization, while in Ref. [18[ the authors use 6tt. 
The results from fitting the tttt channel (dashed curve shown in Fig. [2} are 

m p = 768.6 MeV , g p ^ = 6.59 , T p = 135 MeV , 
while those from fitting the coupled channels (dotted curve shown in Fig. [J]) are 

m p = 769.6 MeV , 5p7r7r = 6.79 , T p = 144 MeV . 
This should be compared with the "exact" results with the chiral unitary approach with two channels 

m p = 767.9 MeV , g p ^ = 6.69 , T p = 139 MeV . 



(16) 

(17) 
(18) 
(19) 
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It is interesting to show explicitly how the levels appear in this case. As shown in Eq. (|13p. the energies are 
obtained for one channel by solving the equation {Vu(E)^ 1 — G U (£,I)) = 0. In Fig. Owe show Vn(E)^ 1 and 
Gn(E,L) separately and the points where the two curves cut each other provide the eigenenergies of the box. We 
observe that G is a small quantity compared with other hadron scattering cases, since it raises very fast close to the 
poles which are given by the free energies of the box induced by the boundary condition on the sides of the box. 
The 7T7T potential is large, and even has a pole, the CDD pole at the bare mass of the p. Consequently Vn(E)^ 1 is 
also a small quantity, passing through zero around the bare p mass. The combination of these two facts leads to an 
eigenenergy of the box with interaction in the vicinity of the bare p mass, which defines the first level. For all the 
other levels the eigenenergies are very close to the free energies and this tells us that they carry small information on 
the potential. All these features are a consequence of the peculiar nature of the p which is not dynamically generated 
by the interaction of the two pions and cannot be described as a scattering resonance from an ordinary potential, but 
is a genuine state (basically a qq state) which decays into two pions. Even then the Liischer formalism, or our version 
of it, can be applied to obtain the p properties from the ttit phase shifts, as we are doing here. 

1 

0.5 


-0.5 
-1 

400 800 1200 1600 2000 
E [MeV] 

FIG. 5: The dashed curve is the ttty inverse potential Vu(E)~ 1 calculated using the chiral unitary approach 41]. The solid line 
is G. 



IV. DEPENDENCE ON THE PION MASS 

In this section, we study the mass and decay width of the p meson using different pion masses. There has been pio- 
neering work along these lines using effective theories in (5(| , evaluating the p self-energy with infrared regularization, 
but the results depend on four low energy constants that ultimately are fitted to lattice data. 

We define mj to be the physical ir mass, and now is a free parameter. It changes from mj to 3mJ in our study. 
When changing m n , other parameters in our previous formulae, especially the parameter / (/„.) used in Eqs. @, can 
also change at the same time. The variation of / w as a function of has been calculated using the IAM in [571 [58j 
and it compares favorably with lattice QCD calculations [59 6lJ. Using the data shown in Fig. 1 of [57j . and after 
fixing / ff (m^) = 87.4 MeV, the value used in our previous calculation, we find the following fit: 

Urn*) = 84.4McV + 3MeV x . (20) 

m% 

The coupling Gy with the vector meson dominance assumption of [H, [H[ is related to f„ [Hj], as Gy = f-n/y/2. We 
still take the value Gy = 53 MeV of [41] but assume it to be proportional to f n as a function of m w . On the other 
hand, the bare p mass, M p in Eqs. offers the link of the theory to a genuine component of the p meson, not 
related to the pion cloud, and we assume it to be m n independent. 

Having all the input, we can use the same procedure to calculate the energy levels using different values. We 
show the energy levels using m T = 1.5 (left), m n — 2.0 mj (middle) and m n = 2.5 mj (right) in Fig. [6] We note 
that the s-coordinate is expressed in units of m^ 1 , not (raj) -1 . 
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We compare our results with the Lattice results of Ref. [l 



Our Input and Result 


Input and Result in Ref. [181 


m w = 276McV, L = 2.75m- 1 = 1.96fm 


m n = 266MeV, L = 1.98fm 


E 1 = 815McV, E 2 = 1390MeV 


Ei = 812MeV, E 2 = 1431MeV 


8 1 = 135.9°, S 2 = 175.8° 


<5i = 130.56°, S 2 = 146.03° 



We find that the first (lowest) energy level and the extracted phase shift from the two approaches are very similar, 
while those from the second energy level have small differences. 

180 
135 
90 
45 


500 700 900 1100 

E [MeV] 

FIG. 7: The nn phase shifts with different pion masses. The solid curve is obtained using physical pion mass, while the 
short-dashed, middle-dashed, long-dashed curves correspond to = 1.5 m„ = 2.0 m% and m„ = 2.5 mj, respectively. 




Following our previous procedure, we can use these energy levels to obtain the phase shifts. Again we only consider 
the first (lowest) energy level. The results are shown in Fig. [7] The solid curve is obtained using the physical pion 
mass, while the short-dashed, middle-dashed, long-dashed curves correspond to = 1.5 mj, m w = 2.0 mj and 
TOtt = 2.5 mj, respectively. 



TABLE I: The p mass, coupling and decay width obtained according to Eq. 1(160 ■ 



ma (mj) 


1.0 1.25 1.5 1.75 2.0 2.25 2.5 2.75 


m p (MeV) 


769.6 772.8 775.7 780.0 785.0 791.3 798.6 806.6 




6.79 6.56 6.32 6.10 5.92 5.62 5.46 5.25 


r (MeV) 


144 119 93.2 69.9 49.1 29.6 15.1 4.26 
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1.0 1.5 2.0 2.5 3.0 1.0 1.5 2.0 2.5 3.0 1.0 1.5 2.0 2.5 3.0 

m„ [m°] m„ [m„°] m, [m°] 

FIG. 8: The p mass (left), the coupling constant (/ p7 r,r (middle) and the decay width T p (right) as functions of m^. 



Using the phase shifts shown in Fig. [71 we can obtain the p mass (left), the coupling constant g p7r7T (middle) and the 
decay width T p (right), which are shown in Fig. [5] as functions of m^. We also show the numbers in Table Q] Again 
we compare our results with the Lattice results in Ref. [Hj]. They use m T = 266 MeV, and the obtained results are 
to p = 792 ± 7 ± 8 MeV and — 5.13 ± 0.20. After correcting for the factor the obtained results in [l8[ 
are m p ~ 792 ± 7 ± 8 MeV and g pffw = 5.93 ± 0.24 in our normalization. The obtained p mass 792 ± 7 ± 8 MeV is 
in agreement, within errors, with our result w p (m T = 274MeV) = 785.0 MeV {n = 2 in Table [J), while the obtained 
coupling constant g^n = 5.93 ± 0.24 is practically identical to our result g^^m^ = 274MeV) = 5.92. 



V. COMPARISON WITH OF OUR RESULTS WITH THE STANDARD LUSCHER'S APPROACH 

In the Lattice QCD calculations, the Liischer's approach is used, which is a bit different from ours as discussed 
in (24l |. To make our analysis complete, especially for the comparison of our results with the Lattice results, in this 
section we show what one would get from the energy levels obtained here using the standard Liischer's approach. It 
was shown in [24| that the exact form of the two meson propagator I(q) of Eq. (fTTj) . where q is the integration or 
summation variable, could be expressed in the following form 

1 wi +u 2 _ _L 1 1 1 , 21 ^ 

2wiw 2 E 2 — (wi + U2) 2 + it 2Ep 2 — cp+ie 2u>iu> 2 wi + W2 + E 

11 11 

4wiW2 ^1 — ^2 — E ALO1U2 UJ2 — U\ — E ' 

where p = X 1 ^ 2 (E 2 , m 2 , m 2 )/2E and X(x, y, z) — (x — y — z) 2 — 4yz. In the standard Liischer's approach only the 
first term of the right hand side of this equation is kept. This guarantees that the imaginary part of the propagator is 
kept exactly, but real terms (which are actually exponentially suppressed in L) are neglected. We want to show here 
the effect of neglecting these terms in the analysis that we have done. 

We obtain Liischer's approach using our formalism, changing I(q) by the first term of Eq. (f2"Tj) . We use as input the 
energy levels which we have calculated in Sec. IIIII using our G- function (Eqs. ([5]) and (flT))) ), and then obtain the phase 
shifts as discussed before but using our expression for I(q) or the approximate one that leads to Liischer's approach. 
Since we have two sets of energy levels, for the single irir channel and for the coupled channels, we take these two sets. 
The results are quite similar for the single and two channel cases and we show the results in Fig. |H1 The left and right 
figures are for single channel and coupled channels "synthetic data", respectively, extracted with just the irir channel 
analysis. We find that in the low energy region the results using the two different approaches are quite similar, while 
the differences become larger as the energy gets bigger. The new results with the Liischer's approach using the single 
7T7T "synthetic data" are 

m p = 786.6 MeV , 5p7r7r = 7.66 , V = 189 MeV , (22) 

while those obtained from the coupled channels "synthetic data" are 

m p = 789.5 MeV , g p7rK = 7.94 , T = 204 MeV . (23) 

We can see that there are quite some differences from using one or the other approach. One should however realize 
that the differences come because the phase shifts for energies around 900-1000 MeV are obtained in our case from 
data with L rather small, of the order of m^ 1 or smaller. In the work of [18| small sizes are avoided by getting data 
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400 500 600 700 800 900 1000 1100 400 500 600 700 800 900 1000 1100 

E [MeV] E [MeV] 

FIG. 9: The comparison between our result (solid curves) with the standard Luscher's approach (dashed curves), respectively. 
The curves at the left and right hand side are extracted from the first energy level using one-channel and two-channel approach, 
respectively. 



using instead the tttt system in moving frames for which the formalism is also available in [65W68I ] , and in [39l . [69j 
using the formalism of 24]. This avoids having to go to small values of L to get higher energies. Conversely, what 
we are telling here is that one could use smaller lattice sizes, with the consequent economy in the lattice calculations, 
but then the analysis would have to be done with our formalism, since the standard Luscher's approach would not 
provide sufficient accuracy. 



VI. THE INVERSE PROCESS OF GETTING PHASE SHIFTS FROM LATTICE DATA FROM TWO 

LEVELS. 

In this section we study the inverse process of getting phase shifts from Lattice Data [24l l33l. l38l l39l. [551 ] using two 
levels and a parametrized potential. The method was shown to be rather efficient in [24l l33l. l38l I39U55T] (see also [38[ 
for other possible parameterizations which induce a bigger dispersion on the obtained results). The energy levels in 
Fig. [2] and Fig. [3] obtained in Sec IIIII contain more information than the one that one obtains using only one level, as 
we have done so far. In this section we assume that they are "Lattice" inputs (synthetic data) and analyze them. We 
shall try to obtain the F-matrices as well as the phase shifts from these "Lattice" data. The procedure followed here 
also serves to get an idea of the uncertainties of the phase shifts for a given accuracy of the "lattice data". 

Since we found that in the region of interest using just the 7T7t channel provides phase shifts with a high accuracy, 
we do the one channel analysis and obtain the I/-matrix. We use the following function which accounts for a CDD 
pole 50]) to do the one-channel fitting: 

V^(s)=cipl(l + -^). (24) 

where c\, C2 and C3 are three free parameters which we shall fit with the "Lattice" data shown in Fig. [5J We note 
that these data are calculated using the coupled channels, but, as mentioned, we find that the single 7T7t channel can 
already give a good enough fitting. We should mention that with this potential and the formalism of the paper we 
can accommodate forms of the T-matrix which are more general than just a Breit-Wigncr, which is usually assumed 
in many Lattice calculations. Our procedure follows Refs. [24|, [Hj]. We take five energies from the first level and five 
more from the second one, and we associate to them an error of 10 MeV. The following set of parameters leads to a 
minimum Xmin = 0.15: 

ci = 9.2 x 10- 5 McV~ 2 , c 2 = 0.73 , c 3 = 843MeV. (25) 

In Fig. [TU] we show our fitting results calculated from all the possible sets of parameters having x 2 < Xmin + 1- We 
find that when the energy becomes larger, the error bars become larger at the same time, but the fitting is quite good. 
Nevertheless, we observe that the percentage in the errors of the induced phase shifts is larger than the one in the 
value of the energies of the box. 

In the former analysis we have used the same regularization that was used in the direct problem of getting phase 
shifts from the chiral unitary approach. However, in a blind analysis one does not know how the G-function should 
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be regularized. As mentioned in [2J| the results of this inverse analysis do not depend on which cut off, or subtraction 
constant one uses in the analysis, as far as one uses the same ones to induce V from the lattice data and then later 
on to get the phase shifts in the infinite volume from Eq. ([!]). For this purpose we change the constant -2 in Eq. (jHJ 
by -1 and -3 and repeat the procedure. The results are shown in Fig. [HI and [T2l As one can see, the results for the 
phase shifts are essentially the same, confirming once more the findings of [24j . 




FIG. 10: Generated data points (10 MeV error, 10 points for each set). The solid curves are the orginal results, and the dashed 
curves (shaddy band) are the fitted results. Fits that fulfill the Xmm + 1 criterion are also shown in all figures (bands). The 
G-function of Eq. © is used here. The discrete points in the figure to the right are the results of the direct determination from 
each "data" points on the figure to the left using Eqs. (|14til5p (see text). 




FIG. 11: Same as in Fig 1101 Here we replace the constant -2 of Eq. Q by — 1. 



In Figs . \TU\ im and [T2l we see the band of values for the phase shifts and some discrete values with error bars. The 
band of values has been obtained with a fit to the synthetic data with a potential as explained before (fit method). 
The discrete points correspond to the extraction via Eqs. (| 1 411 1 5 [) of the phase shifts from each individual synthetic 
data of the five energies of the lower level shown on the left of the figures, with their corresponding error (direct 
method). We can see that the errors in the phase shifts are large at small energies but they become smaller as the 
energy increases. Yet, the overall fit to the whole set of the 10 data from two levels provide a more accurate global 
determination of the phase shifts than the direct method in the whole energy range. 

VII. CONCLUSION 

We have made a study of the irir interaction in p-wave in a finite box using the chiral unitary approach which is 
very successful to provide phase shifts and the pion form factor. With this procedure we obtain discrete energies in 
the box as a function of L and of which compare favorably with present lattice results. The second part consist of 
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FIG. 12: Same as in Fig 1101 Here we replace the constant -2 of Eq. © by —3. 



getting the tttt phase shifts in infinite volume from these finite box energies for what we use our equivalent procedure 
to Liischer's approach developed in [iH, which improves it for small lattice sizes. We use two methods, the first one 
(direct method), valid for only one channel, produces the phase shifts without having to use an intermediate potential. 
The second one uses a parametrized potential and makes a fit to data. This second procedure takes advantage of the 
correlation among data and provides a global fit that allows one to get phase shifts for all energies and not only for 
the energies which are eigenvalues of the box, as in the first method, or in Liischer's approach. 

We also estimated errors induced in the phase shifts from the errors in the energies of the box. We could see that by 
taking a set of about 10 energies from two levels one could induce good phase shifts with reasonable errors, although 
we noted that the percentage of the errors in the phase shifts was larger than the one in the energies, Typically, errors 
of 10 MeV in the eigenenergies of the box resulted in large errors for small energies using the direct method but smaller 
errors at energies around 800 MeV. The fit method provides a better global agreement for all energies and the errors 
in the phase shifts, induced by the 10 MeV uncertainties in the eigenenergies of the box, are of the order of 6%. 

So far we have used the centroid of the synthetic data centered at the exact curve. In actual lattice data there 
will also be a dispersion of the centroids. In [24( this was taken into account and a good statistical study was done. 
The conclusion, which we can accept also in the present case, is that taking into account a dispersion of the centroids 
of about 5 MeV, produces extra uncertainties in the results of maximum 50% of those produced by the 10 MeV 
uncertainty of the data. Hence we end with a round 10% errors in our case. 

We also showed that if one uses periodic boundary conditions in the rest frame of the irir system one has to use 
small values of the box size where the standard Liischer approach produces too large errors. While this problem has 
been avoided in some lattice works by obtaining data in a 7T7t moving frame, we show that using our procedure one 
can go to smaller sizes with a considerable saving of the computer time in actual QCD lattice calculations. These 
findings can serve as a guideline for future studies of the p meson in these calculations, which would also allow to use 
pion masses closer to the physical one. 

Acknowledgments 

We would like to thank B. X. Sun for help in the calculations at the early stage of the problem. A careful reading 
and comments from U.-G. Meifiner, A. Martinez Torres, M. Doring and L. Roca are much appreciated. This work is 
partly supported by DGICYT contract FIS2011-28853-C02-01, the Generalitat Valenciana in the program Prometeo, 
2009/090, and the EU Integrated Infrastructure Initiative Hadron Physics 3 Project under Grant Agreement no. 
283286, the National Natural Science Foundation of China under Grant No. 11147140, and the Scientific Research 
Foundation for the Returned Overseas Chinese Scholars, State Education Ministry. 



[1] Y. Nakahara, M. Asakawa, T. Hatsuda, Phys. Rev. D60 (1999) 091503. 

K. Sasaki, S. Sasaki and T. Hatsuda, Phys. Lett. B 623 (2005) 208. 
[2] N. Mathur, A. Alexandra, Y. Chen et al, Phys. Rev. D76 (2007) 114505. 



12 



[3] S. Basak, R. G. Edwards, G. T. Fleming et al, Phys. Rev. D76 (2007) 074504. 
[4] J. Bulava, R. G. Edwards, E. Engelson et al, Phys. Rev. D82 (2010) 014507. 
[5] C. Morningstar, A. Bell, J. Bulava et al, AIP Conf. Proc. 1257 (2010) 779. 
[6] J. Foley, J. Bulava, K. J. Juge et al, AIP Conf. Proc. 1257 (2010) 789. 
[7] M. G. Alford and R. L. Jaffe, Nucl. Phys. B 578 (2000) 367. 

[8] T. Kunihiro, S. Muroya, A. Nakamura, C. Nonaka, M. Sekiguchi and H. Wada [SCALAR Collaboration], Phys. Rev. D 

70 (2004) 034504 

[9] F. Okiharu et al, |arXiv:hep-ph/0507i87] H. Suganuma, K. Tsumura, N. Ishii and F. Okiharu, PoS LAT2005 (2006) 070; 

Prog. Theor. Phys. Suppl. 168 (2007) 168. 
[10] C. McNeile and C. Michael [UKQCD Collaboration], Phys. Rev. D 74 (2006) 014508; A. Hart, C. McNeile, C. Michael 

and J. Pickavance [UKQCD Collaboration], Phys. Rev. D 74 (2006) 114504. 
[11] H. Wada, T. Kunihiro, S. Muroya, A. Nakamura, C. Nonaka and M. Sekiguchi, Phys. Lett. B 652 (2007) 250. 
[12] S. Prelovsek, C. Dawson, T. Izubuchi, K. Orginos and A. Soni, Phys. Rev. D 70 (2004) 094503; S. Prelovsek, T. Draper, 

C. B. Lang, M. Limmer, K. F. Liu, N. Mathur and D. Mohler, Conf. Proc. C 0908171 (2009) 508; Phys. Rev. D 82 (2010) 

094507. 

[13] H. -W. Lin et al [ Hadron Spectrum Collaboration ], Phys. Rev. D79, 034502 (2009). 

[14] C. Gattringer, C. Hagen, C. B. Lang, M. Limmer, D. Mohler, A. Schafer, Phys. Rev. D79, 054501 (2009). 

[15] G. P. Engel et al [ BGR [Bern-Graz-Regensburg] Collaboration ], Phys. Rev. D82, 034505 (2010). 

[16] M. S. Mahbub, W. Kamleh, D. B. Leinweber, A. O Cais, A. G. Williams, Phys. Lett. B693, 351-357 (2010). 

[17] R. G. Edwards, J. J. Dudek, D. G. Richards, S. J. Wallace, Phys. Rev. D84, 074508 (2011). 

[18] C. B. Lang, D. Mohler, S. Prelovsek and M. Vidmar, Phys. Rev. D 84, 054503 (2011). 

[19] S. Prelovsek, C. B. Lang, D. Mohler and M. Vidmar, PoS LATTICE2011, 137 (2011). 

[20] M. Liischer, Commun. Math. Phys. 105 (1986) 153. 

[21] M. Liischer, Nucl. Phys. B 354 (1991) 531. 

[22] V. Bernard, U. -G. Meissner, A. Rusetsky, Nucl. Phys. B788, 1-20 (2008). 

[23] V. Bernard, M. Lage, U. -G. Meissner, A. Rusetsky, JHEP 0808, 024 (2008). 

[24] M. Doring, U. -G. Meissner, E. Oset and A. Rusetsky, Eur. Phys. J. A 47, 139 (2011). 

[25] X. Feng et al [ETM Collaboration], PoS LAT 2010, 104 (2010). 

[26] S. Aoki et al. [CP-PACS Collaboration], Phys. Rev. D 76, 094506 (2007). 

[27] M. Gockeler et al [QCDSF Collaboration], PoS LATTICE 2008, 136 (2008). 

[28] S. Aoki et al [PACS-CS Collaboration], PoS LATTICE 2010, 108 (2010). 

[29] X. Feng, K. Jansen and D. B. Renner, Phys. Rev. D 83, 094505 (2011). 

[30] J. Frison et al [Budapest-Marseille- Wuppertal Collaboration], PoS LATTICE 2010, 139 (2010). 

[31] C. McNeile et al [UKQCD Collaboration], Phys. Lett. B 556, 177 (2003). 

[32] M. Doring, J. Haidenbauer, U. -G. Meissner, A. Rusetsky, Eur. Phys. J. A 47, 163 (2011). 

[33] A. Martinez Torres, L. R. Dai, C. Koren, D. Jido and E. Oset, Phys. Rev. D 85, 014027 (2012). 

[34] E. E. Kolomeitsev, M. F. M. Lutz, Phys. Lett. B582, 39-48 (2004). 

[35] J. Hofmann, M. F. M. Lutz, Nucl. Phys. A733, 142-152 (2004). 

[36] F. -K. Guo, P. -N. Shen, H. -C. Chiang, R. -G. Ping, B. -S. Zou, Phys. Lett. B641, 278-285 (2006). 

[37] D. Gamermann, E. Oset, D. Strottman, M. J. Vicente Vacas, Phys. Rev. D76, 074016 (2007). 

[38] M. Doring, U. G. Me issner, JH EP 1201, 009 (2012). 

[39] L. Roca and E. Oset. larXiv:1201.0438l [hep-lat], 

[40] J. A. Oiler and E. Oset, Phys. Rev. D 60, 074023 (1999). 

[41] J. A. Oiler, E. Oset and J. E. Palomar, Phys. Rev. D 63, 114009 (2001). 

[42] J. Gasser and H. Leutwyler, Nucl. Phys. B 250, 465 (1985). 

[43] G. Ecker, J. Gasser, A. Pich and E. de Rafael, Nucl. Phys. B 321, 311 (1989). 

[44] A. Dobado and J. R. Pelaez, Phys. Rev. D 56, 3057 (1997). 

[45] J. A. Oiler, E. Oset a nd J. R. Pelaez, P hys. Rev. D 59, 074001 (1999) [Erratum-ibid. D 60, 099906 (1999)] [Erratum- ibid. 

D 75, 099903 (2007)] hep-ph/9804209] . 
[46] B. Borasoy and R. Nissler, Eur. Phys. J. A 26, 383 (2005). 
[47] J. R. Pelaez, Phys. Rev. Lett. 92, 102001 (2004). 
[48] J. R. Pelaez and G. Rios, Phys. Rev. Lett. 97, 242002 (2006). 
[49] J. Nieves and E. Ruiz Arriola, Phys. Rev. D 80, 045023 (2009). 
[50] L. Castillejo, R. H. Dalitz and F. J. Dyson, Phys. Rev. 101, 453 (1956). 

[51] J. A. Oiler and E. Oset, Nucl. Phys. A 620, 438 (1997) [Erratum-ibid. A 652, 407 (1999)] [hep-ph/9702314] . 

[52] J. A. Oiler and U. G. Meissner, Phys. Lett. B 500, 263 (2001). 

[53] P. Estabrooks and A. D. Martin, Nucl. Phys. B 79, 301 (1974). 

[54] J. Gasser and H. Leut wyler, Nucl. Phys . B 250, 517 (1985). 

[55] J. J. Xie and E. Oset. larXiv:1201.0149l [hep-ph], 

[56] P. C. Bruns and U. G. Meissner, Eur. Phys. J. C 40, 97 (2005). 

[57] J. R. Pelaez and G. Rios, Phys. Rev. D 82, 114002 (2010). 

[58] C. Hanhart, J. R. Pelaez and G. Rios, Phys. Rev. Lett. 100, 152001 (2008). 

[59] P. .Boucaud et al [ETM Collaboration], Phys. Lett. B 650, 304 (2007). 

[60] S. R. Beane, T. C. Luu, K. Orginos, A. Parreno, M. J. Savage, A. Torok and A. Walker-Loud, Phys. Rev. D 77, 014505 



13 



(2008). 

[61] J. Noaki, S. Aoki, T. W. Chiu, H. Fukaya, S. Hashimoto, T. H. Hsieh, T. Kaneko and H. Matsufuru et ai, PoS LATTICE 
2008, 107 (2008). 

[62] J.J. Sakurai, Currents and mesons (University of Chicago Press, Chicago II 1969). 
[63] M. Bando, T. Kugo and K. Yamawaki, Phys. Rept. 164, 217 (1988). 

[64] G. Ecker, J. Gasser, H. Leutwyler, A. Pich and E. de Rafael, Phys. Lett. B 223, 425 (1989). 

[65] K. Rummukainen, S. A. Gottlieb, Nucl. Phys. B450, 397-436 (1995). 

[66] C. h. Kim, C. T. Sachrajda, S. R. Sharpe, Nucl. Phys. B727, 218-243 (2005). 

[67] Z. Davoudi and M. J. Savage, Phys. Rev. D 84, 114502 (2011). 

[68] Z. Fu, Phys. Rev. D 85, 014506 (2012). 

[69] M. Doring, U.G. Meissner, E. Oset and A. Rusetsky, to be submitted. 



